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EXPERIMENTAL  UNCERTAINTY  IN  THIRD-ORDER  SUSCEPTIBILITY 
AND  MOLECULAR  SECOND  HYPERPOLARIZABILITY  MEASUREMENTS 
OBTAINED  WITH  DEGENERATE  FOUR  WAVE  MIXING 

INTRODUCTION 

Third-order  susceptibility  and  molecular  second  hyperpolarizability 
measurements  are  currently  performed  with  degenerate  four-wave  mixing  (DFWM) . 
Experimental  values  obtained  by  this  technique  have  been  reported  with 
estimated  confidence  intervals*  This  report  describes  an  uncertainty  analysis 
that  can  be  used  to  1)  better  define  experimental  uncertainties,  2)  add  some 
insight  into  the  degree  to  which  individual  experimental  parameters  affect 
overall  system  accuracy,  and  3)  validate  comparisons  of  experimental 
(y)  values  for  various  compounds  to  confirm  the  significance  of  any  observed 
variations  in  nonlinear  response,  as  in  the  study  of  trends* 

The  DFWM  system  described  here  utilizes  a  backward-beam  geometry  where 
two  strong  collinear  pump  beams,  propagating  in  opposite  directions,  interact 
at  the  sample  with  a  third  less-intense  probe  beam  to  create  a  fourth 
nonlinear  optical  phase  conjugate  beam.  The  intensity  of  the  conjugate  signal 
is  proportional  to  the  square  of  the  x^^^  of  the  sample  under  study.  The 
signal  is  measured  with  a  silicon  photodiode  and  oscilloscope*  It  is  then 
compared  to  a  similar  signal  obtained  from  a  carbon  disulfide  standard  under 
identical  conditions  using  the  expression 


where  n  is  the  index  of  refraction,  a  is  the  linear  absorption  coefficient,  d 
is  the  length  of  the  sample  and  I  is  the  measured  signal  of  the  phase 
conjugate  beam*  This  approach  eliminates  the  need  of  obtaining  absolute 
irradiance  measurements  for  each  of  the  additional  three  beams  used  in  the 
DFWM  setup*'’ 

The  relationship  between  the  molecular  second  hyperpolarizability  (y) 
and  the  third-order  macroscopic  susceptibility  of  the  solute  is  given  by  the 
expression 


(2) 


where  N  is  the  number  density  of  molecules  per  rtiL  and 


L  = 


n^+2 

3  ■ 


(3) 


It  is  generally  assumed  that  the  conjugate  signal  is  due  to  a 
combination  of  third-order  nonlinearities  in  both  the  solvent  and  solute  and 
that  for  a  solution  of  noninteracting  particles  the  relationship  is  addrtrve: 


solvtion  ^\^solv0nty solvent* ^soluce^ soJute] ' 


(4) 


The  X^^^solution  obtained  experimentally  with  Eg.  1.  For  nonsaturated 
solutions 


solute 


=  A£ 

M 


(5) 


where  A  is  Avogadro's  number,  M  is  the  molecular  weight  of 
is  the  concentration.  Making  appropriate  substitutions  of 
Eg.  4,  the  following  expression  is  obtained: 

X^ solution  ~  X^ solvent  *  y solute 

By  evaluating  Eg.  1  for  a  series  of  dilute  solutions  and  plotting  the 
results  versus  respective  concentration,  a  straight  line  is  obtained.  A 
least— sguares  approximation  is  then  applied  to  the  data  from  which  one  obtains 

X<^>solvent  y- intercept  and  Ygolute  slope.  The  X^^^olute 

obtained  from  Eg.  2. 


the  solute,  and  C 
Egs.  2  and  5  into 

(6) 
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DERIVATION  OF  UNCERTAINTY  EQUATIONS 


If  it  is  assumed  that  the  uncertainties  in  the 

independent  variables  (X^,  X2,...,X„)  are  all  given  with  the  same  probability, 
then  the  uncertainty  (1(7^)  in  the  dependent  variable  (R)  is  given  by  the 
expression 


(7) 


By  inserting  Eg.  1  into  Eg.  7  a  probable  estimate  of  the  uncertainty  in 
experimentally  determined  X^^^sanple  obtained: 


where 


s  _  (  uL  ^ah/z 


(9) 


ixfL.  _ 

51*  \  j  I  ^s] 


(tL  ^qlL/2 


(10) 


5X»>s  _  fia)  «!.  e.^/2 

5J«  2  (j^)3/2  (^d*j  (i-e-“^) 


(11) 


y(3)  ^  2^*"^  / 


aL  ^aL/2 

■■■■  "  # 


(l-e-«i) 


(12) 


-  v(3)  f-^1 


clL  ^«L/2 


(13) 


3 


=  v3)  fifV  ffsf  (± 


^«L/2 

(l-e-«^)  ' 


<14) 


aj^u) 


dd. 


S  _  „(3> 

X  R 


®(tF(^) 


aL 


il-e-^) 


j«X./2 


(15) 


and 


3^(3)^ 

da 


(16) 


Similarly,  for  the  uncertainty  in  (y)»  Eq.  2  is  applied  to  Eg.  7  to 
obtain: 


ffy  = 


1/2 


(17) 


where 


ay  ^  _iL 
ax‘^’ 


(18) 


&Y  - 

dM  l*A 


(19) 


ay  =  -4^<3>jf 

ai,  1,5^ 


(20) 


and 


=  0  (assumed)  . 
dA 


(21) 


In  order  to  evaluate  Eq.  8^  absolute  error  values  are  recjuired  for 
Jg,  n^r  n^,  d^,  dg,  and  a.  These  can  all  be  determined  directly  through 


4 


experimental  considerations  or  can  be  given  assumed  values.  For  Eq.  17  error 
values  for  K,  (solute),  and  L  must  be  obtained.  Additional  calculations 

must  be  performed  to  obtain  estimates  for  the  uncertainty  in  (solute), 

and  L.  For  W^,  Eqs  3  and  7  are  combined; 


= 


(22) 


where 


dL  _ 
dn,  ■  3 


(23) 


For  the  uncertainty  in  x^^^/  ^  first  written  in  the  form; 


(3>  «  v(3) 

s 


(1000) 


(24) 


The  factor  of  1000  is  included  as  a  dimensional  constant.  Applying  Eq.  7; 


A  Mlt 


dx 


U> 


sit 


V  / 


3% 


(3) 


A  i 


ax 


(3) 


sit 


3x 


(3) 


Slv 


*'=) 


1/2 


(25) 


where 


^X^sit  _  1000 

^x^  ■  C  • 


(26) 


H^sit  ^  -1000 


and 


dX^slv 


=  1000 


y(3) 

A  S 


_  v(3) 

X  slv 


(28) 
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The  value  of  depends  upon  the  method  with  which  the  solutions  are 
prepared*  The  technique  currently  being  used  relies  on  accurate  weight 
measurements  to  determine  the  solution  concentration*  Briefly,  an  initial 
concentration  (C- )  is  prepared  using  a  precision  balance  for  the  solute  and  a 
volumetric  flask  for  the  solvent.  For  each  run  a  portion  of  the  solution  is 
used.  After  each  run  the  mass  of  the  remaining  solution  is  determined 
The  flask  is  again  filled  and  weighed  to  obtain  the  initial  weight  of  the  next 
run  (iw^).  Using  the  following  expression  the  new  concentration  is  obtained 
(see  Appendix  A  for  derivation): 


(29) 


Advantages  of  this  technique  include  ease  of  weight  measurements,  a  reduction 
in  waste  solvent,  and  high  overall  accuracy*  To  obtain  we  evaluate  the 
following  equation; 


(30) 


where 


dCi  m„ ' 


(31) 


^  =  £i 

diUf  m„ 


(32) 


and 


^  =  zSiUk. 

dm„  m/ 


(33) 


In  order  to  evaluate  Eq*  30  a  value  for  must  be  estimated*  Beginning  with 
the  expression 


(34) 
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where  m,-  is  the  initial  weight  and  V,-  is  the  initial  volume,  we  must  evaluate 
the  following  equation: 


(35) 


where 


^  -  A 

dnii  ~  Vi 


(36) 


and 


(37) 


These  expressions  have  been  incorporated  into  a  computer  program 

(Appendix  B)  which  calculates  X^^^olution'  X^’^olute'  solute 

uncertainties  associated  with  each  of  these  values. 
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EXPERIMENTAL-PARAMETER  SENSITIVITY 


In  order  to  evaluate  the  uncertainty  equations  that  were  derived  in  the 
previous  section,  values  are  assumed  for  the  absolute  errors  of  the 
independent  variables.  These  are  qiven  as  percent  errors  in  Table  1. 


Table  1.  Assumed  percent  errors  for  experimental  parameters 


Parameter 

Assumed  Percent  Error 

v(3) 

*  CS2 

0% 

^solution 

5% 

■^CS2 

5% 

^solution 

1% 

”CS2 

0.1% 

^sample 

0.1% 

^CS2 

0.1% 

^solution 

0.08% 

^solute 

0.1% 

®^^'®solution 

1% 

voiumesotution 

1% 

It  is  assumed  that  an  exact  value  for  known.  Uncertainties  in 

intensity  measurements  are  based  on  an  estimate  of  the  readability  of  the 

oscilloscope  used  in  the  evaluations.  The  index  of  refraction  for  CS2  is 

considered  more  accurate  than  that  of  the  sample  because  it  is  assumed  that 

n  .  =  n  I  and  additional  error  is  allowed  here.  Percent  errors  given 

solution  solvent' 

for  the  cell  thicknesses  <?sample  ^CS2  manufacturer  stated  tolerances. 
Mass  and  volume  tolerances  are  based  on  scale  readabilities. 

Absolute  errors  for  solution  concentrations  are  calculated  using  Eq.  30. 
These  will  slightly  increase  with  each  new  solution  that  is  prepared.  For 
example,  for  an  initial  concentration  of  1  mg/mL,  where  0.005  g  of  solute  is 
dissolved  in  5  itiL  of  solvent  with  tolerances  of  ±0.0001  g  and  ±  0.05  mL, 
respectively.  Eg.  35  gives  an  uncertainty  of  0.00002  g/mL  or  2.2%  for  the 
initial  concentration.  This  value  is  then  used  as  input  into  Eq.  30  to  obtain 
the  uncertainty  in  the  following,  more  dilute  solution.  Table  2  shows  an 
example  in  which  each  consecutive  concentration  is  half  the  previous 
concentration.  New  concentrations  (C^)  are  calculated  at  the  end  of  each  run 
using  the  mass  measurements  shown  for  that  run  and  Eq.  29.  Absolute 
uncertainties  are  predicted  with  Eq.  30  and  percent  errors  are  also 

shown  for  each  run. 
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Table  2. 

Uncertainty  in  concentration 

measurements 

for  a  standard 

run. 

“i 

“f 

“n 

Cn 

Wen  » 

Error 

(g/mL) 

(g) 

(g) 

(g) 

(g/mL) 

(g/mL) 

IxlO'-* 

4.6 

2.3 

4.6 

5x10"* 

IxlO"’ 

1% 

5x10’^ 

4.6 

2.3 

4.6 

2.5x10'* 

5x10'^ 

1% 

2.5x10'^ 

4.6 

2.3 

4.6 

1.25x10'* 

4x10'* 

1.6% 

1.25x10'* 

4.6 

2.3 

4.6 

6.25x10'® 

2x10'* 

1.6% 

6.25x10'5 

4.6 

2.3 

4.6 

3.12x10'® 

IxlO'* 

1.6% 

These 

values  are 

used  as  input 

to  program  ERRCHI- 

■3. BAS.  Table  3 

gives 

additional 

input  for  a 

representative 

run 

where  the  of  zinc  tetraphenyl 

tetrabenzporphyrin  in  tetrahydrafuran  was  evaluated. 


Table  3.  ERRCHI~3.BAS  input  data  for  experimental  run  of  zinc  tetraphenyl 


tetrabenzporphyrin  in  tetrahydrafuran. 


Variable 

Input 

^solution 

28000 

^CS2 

150000 

oc 

0.875 

0.00076  g/mL 

y(3) 

^  solvent 

2.4x10’^^  esu 

^X<3)solvent 

9.7x10’'*^  esu 

M 

678 

weight 

0.0038  g 

vol 

5  mL 

The  following  output  was  obtained  for  this  run: 

solution  =  1.45x10“^^ 

Uncertainty  =  5.94x10“^^ 

Percent  Error  =  4% 

of  solute  =  1.88x10’^ 

Uncertainty  =  T.USxlO’^"* 

Percent  Error  =  4% 


(y)  of  solute  =  6.3x10"^^ 

Uncertainty  =  3*64x10'^^ 

Percent  Error  =  6% 

Uncertainty  in  X^^^olution 
y<3)  =  0 

Voltage  reading  for  solution  =  3.6x10'^^ 

Voltage  reading  for  CSj  =  3.6x10*^^ 

Solution  refractive  index  =  2.9x10’^^ 

CSj  refractive  index  =  2.9xl0‘^® 

Solution  cell  thickness  =  1.5x10*^^ 

CSj  cell  thickness  =  1.5x10"^^ 

Scunple  absorption  =  5.4x10*^^ 

Uncertainty  in  concentration  due  to: 

Weight  =  0.0076 
Volume  =  0.0076 

Uncertainty  in  X^^^solute 

-  7.8X10-" 

<“’sol,.o.  -  1.3X10-" 

Concentration  =  2.7x10“^^ 

Uncertainty  in  field  factor  due  to: 

Solution  refractive  index  =  0.014 

Uncertainty  in  (y)  due  to: 

=  2-7xl0-» 

Molecular  weight  =  6.3x10*^ 

Field  factor  (L)  =  2.5x10'^^ 

The  value  for  X^^^soivent  corresponding  uncertainty  was  obtained  by 

first  running  the  program  at  zero  concentration.  Uncertainties  due  to 
individual  experimental  parameters  shown  above  are  obtained  from  the  product 
of  the  parameter  sensitivity  and  its  absolute  uncertainty,  for  instance,  from 
Eg.  7 : 

dR  „  (38) 
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These  values  show  the  relative  degree  to  which  single  variables  contribute  to 
the  overall  uncertainty.  For  excunple/  the  greatest  degree  of  uncertainty  in 
the  experimental  value  for  X^^^solution  oscilloscope  voltage  readings. 

A  similar  output  to  that  shown  above  is  obtained  for  each  concentration 

evaluated  in  the  series.  Each  value  of  X^^^solution'  ^C^^^solute'  will  have 

an  associated  uncertainty  based  on  the  values  of  the  experimental  parameters 
unique  to  that  run.  The  final  values  reported  for  X^^^solution'  ^C^^^solute' 

(y)  are  obtained  via  a  curve  fitting  technique  while  their  respective 
uncertainties  are  obtained  by  taking  the  sample  standard  deviation  of  the 
individual  uncertainties  obtained  for  each  run. 

In  some  instances  the  input  uncertainties  given  in  Table  1  may  vary; 
for  example,  in  a  case  where  the  sample  is  in  the  fom  of  a  thin  film,  the 
thickness  of  which  can  only  be  approximated.  With  this  in  mind  a  series  of 
plots  were  generated  which  show  the  change  in  the  uncertainty  of  X^^^solute 
(y)  resulting  from  an  increase  in  the  uncertainty  of  a  single  variable.  These 
are  shown  in  Appendix  C,  Figs.  C1-C21.  In  each  plot  the  uncertainty  in  a 
single  independent  variable  has  been  increased  from  0%  to  50%  error  while  the 
uncertainty  in  the  remaining  independent  variables  has  been  kept  constant  at 
the  current  values. 

The  series  of  sensitivity  plots  for  both  X^^^solute  summarized 

in  the  following  bar  charts  where  the  normalized  average  slope  of  each 
variable  is  shown.  The  charts  show  the  relative  degree  to  which  a  given 
variable  is  sensitive  to  changes  in  percent  error.  For  example,  the 
uncertainty  in  X^^^solute  most  sensitive  to  uncertainty  in  the  refractive 
indices  of  both  CS2  and  the  solution  while  accuracy  in  determining  the  linear 
absorption  of  the  sample  has  a  minimal  effect. 
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Fig.  1.  Normalized  average  slopes  of  parameter  sensitivity 

curves  for  X^^^solute’ 


Parameters  shown  in  Fig.  1: 

1  —  Sample  absorption 


solution 


2  - 
3  ~  Ir 


4  __  y(3) 

^  X  CS2 


5  —  Sample  thickness 


6  —  CSj  thickness 

7  —  Sample  mass 

8  —  Sample  volume 

5  •"cs2 

10  ■^solution 


Figure  2  shows  that  {y)  is  also  most  sensitive  to  uncertainty  in  the 

refractive  indices  of  CSg  and  sample  solution.  The  increased  dependence  on 

n  ,  is  due  to  its  presence  in  the  local  field  factor  (Equation  3). 
solution 
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Normalized  Slope 


1  2345  6789  10  11 


Experimental  Parameter 


Fig.  2.  Normalized  average  slopes  of  parameter  sensitivity 

curves  for  { y) . 


Parameters  shown  in  Fig.  2: 

1  —  Sample  absorption 

2  “  T 

solution 

^  -^csa 

■  4  __  y(3) 

^  X  CS2 

5  —  Sample  thickness 


6  —  CS2  thickness 


7  —  Molecular  weight 

8  —  Sample  mass 

9  —  Sample  volume 

10~  iir 


11—  n 


^CS2 

solution 
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RESULTS 


A  sampling  of  eight  compounds  evaluated  for  ^y)  reveals  average 

X^^^  and  (y>  uncertainties  of  8%  and  9%,  respectively,  with  values  ranging 
from  6%  to  10%  for  x^^^  and  6%  to  12%  for  (y).  In  general,  the  lowest  percent 
error  is  attributable  to  the  highest  voltage  signal.  Higher  x^^^  and  (y) 
values  have  higher  associated  phase-conjugate  voltage  readings  and,  therefore, 
lower  percent  uncertainties.  For  example,  a  material  with  a  x^^^  area 

of  10’^  will  have  a  percent  uncertainty  of  about  7%  whereas  a  x^^^ 
have  an  uncertainty  as  high  as  40%. 

Table  4  shows  the  results  obtained  for  zinc  tetraphenyl 
tetrabenzporphyrin  in  tetrahydrafuran.  It  can  be  seen  from  these  data  that 
while  the  absolute  uncertainty  varies  with  concentration,  the  percent 
uncertainty  remains  relatively  constant. 


Table  4.  Experimental  results  for  zinc  tetraphenyl 
tetrabenzporphyrin  in  tetrahydrafuran. 


c 

(g/mL) 

^  solution 
(xlO"^^  esu) 

^  solute 
(xlO*^^  esu) 

<Y> 

(xlO'^’  esu) 

0.00216 

43.7 

±1.78 

(4%) 

20.1  ±0.82 

(4%) 

6.66 

±0.38 

(6%) 

0.00076 

14.5 

±0.59 

(4%) 

18.8  ±0.78 

(4%) 

6.26 

±0.36 

(6%) 

0.00038 

9.57 

±0.39 

(4%) 

24.5  ±1.04 

(4%) 

8.12 

±0.48 

(6%) 

0.00019 

6.34 

±0.25 

(4%) 

32.1  ±1.44 

(4%) 

10.6 

±0.64 

(6%) 

0.00010 

2.31 

±0.09 

(4%) 

21.2  ±1.14 

(5%) 

7.03 

±0.48 

(7%) 

0 

0.24 

±0.01 

(4%) 

— 

— 

“ 

Figure  3  shows  the  curve  fit  data  from  Table  4  where  X^^^golute  versus 
concentration  is  shown.  The  straight  line  is  obtained  by  plotting  Eg.  24  with 
X^^^solvent  ~  2.4x10’^^  esu,  X^^^solute  “  2.0x10“^  esu,  and  the  concentrations  shown 
in  Table  4.  Specific  values  for  X^^^solution  compared  in  Table  5 
to  experimental  values  obtained  from  Table  4.  Percent  errors  given  in  Table  5 
are  calculated  with  the  assumption  that  the  curve-fit  values  are  the  more 
correct.  It  is  generally  expected  that  experimental  deviations  from  the 
curve-fit  line  will  be  within  one  or  two  standard  deviations  as  given  in  Table 
4,  or  at  most  about  10%.  Absolute  errors  shown  in  Table  5,  in  two  cases,  are 
outside  predicted  uncertainty.  It  is  currently  believed  that  these  relatively 
large  percent  errors  are  due  to  fluctuations  in  laser  power;  a  variable  which 
has  not  been  considered  in  the  uncertainty  analysis.  Variations,  however, 
appear  to  be  random  and  are  expected  to  average  out  in  the  curve  fitting 
process. 
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Fig.  3.  Experimental  data  and  curve-fit  line  versus  concentration  for 
zinc  tetraphenyl  tetrabenzporphyrin  in  tetrahydrafuran. 


Table  5.  Comparison  of  experimental  and  curve-fit  data  for 


zinc,  tetraphenyl  tetrabenzporphyrin  in  tetrahydrafuran. 


c 

(g/mL) 

Experimental 

y<3) 

*  solution 
(xlO*^^  esu) 

Curve-Fit 

y(3) 

^  solution 
(xlO"^^  esu) 

Percent 

Error 

0.00216 

43.7 

43 

2% 

0.00076 

14.5 

15 

4% 

0.00038 

9.57 

7.8 

23% 

0.00019 

6.34 

4.0 

57% 

0.00010 

2.31 

2.2 

5% 

0 

0.24 

0.24 

0% 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Uncertainty  expressions  have  been  derived  and  used  to  evaluate  a  DFWM 
scheme  used  in  the  determination  of  third-order  susceptibility  and  molecular 
second  hyperpolarizability  of  nonlinear  optical  materials. 

The  results  of  a  representative  experimental  run  are  presented  showing 
that,  for  samples  in  solution,  the  greatest  uncertainty  is  currently 
introduced  by  phase-conjugate  voltage  readings.  As  the  strength  of  the  phase- 
conjugate  signal  decreases,  the  overall  experimental  uncertainty  increases. 
This  relationship  can  occur  as  a  result  of  a  low  laser  output  intensity  or  by 
the  nonlinear  material  having  a  low  phase-conjugate  signal  strength.  Maximum 
experimental  accuracy  can  be  assured  with  proper  laser  maintenance  and  system 
alignment  to  optimize  laser  energy  output  and  phase-conjugate  signal  strength. 

Parameter  sensitivity  curves  reveal  that  the  indices  of  refraction  of 
both  the  sample  in  solution  and  reference  solution  are  the  most  critical 
parameters  with  respect  to  overall  experimental  uncertainty.  Solution 
refractive  indices  are  currently  assumed  equal  to  those  of  the  solution 
solvent;  the  validity  of  this  assumption  should  be  verified  in  light  of  the 
results  of  this  study. 

Sensitivity  plots  contained  in  this  report  can  be  used  as  a  general 
guide  in  estimating  the  impact  of  changes  in  tolerances  of  experimental 
variables.  Accompanying  software  is  now  routinely  used  in  DFWM  experiments  in 
determining  and  (y)  and  their  associated  uncertainties. 

Actual  experimental  deviations  from  a  least-squares  curve  were  shown  to 
be  within  predicted  uncertainties  with  the  exception  of  some  outliers 
attributable  to  variation  in  laser  output  power,  which  was  not  considered  in 
the  analysis. 

.  Concerning  the  investigation  of  trends:  results  show  that  changes  of 
±10%  percent  or  greater  in  x^^^  aud  (y)  can  be  considered  statistically 
significant.  This  confidence  interval  can  be  decreased  by  increasing  the 
number  of  points  (concentrations)  used  in  the  least-squares  approximation. 

The  method  of  determining  solution  concentration  outlined  here  has 
proven  to  be  of  acceptable  accuracy  (±1.6%)  while  also  minimizing  the 
production  of  hazardous  waste  and  the  time  required  to  complete  an 
experimental  run. 

It  is  believed  that  this  study  has  given  significant  insight  into  the 
experimental  uncertainties  inherent  to  the  DFWM  technique  when  used  in  the 
evaluation  of  third-order  susceptibility  and  molecular  second 
hyperpolarizability  of  nonlinear  optical  materials. 
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APPENDIX  A.  Derivation  of  Mass/Concentration  Expression 


A  technique  has  been  developed  to  determine  dilute  solution 
concentrations  from  accurate  weight  measurements.  The  derivation  is  as 
follows: 


An  initial  predetermined  concentration  C-  is  prepared  of  weight  and  density 
p.  Some  solution  is  removed  for  measurement  and  the  sample  is  weighed  again. 
This  final  weight  is  jn^.  To  determine  the  final  volume  we  can  write: 


Vf  = 


9  ‘ 


(Al) 


The  mass  of  the  solute  contained  in  can  be  determined  by  the  expression 


^Bolvte  “  Tf* 


(A2) 


Some  solvent  is  added  to  dilute  the  solution  for  the  next  measurement.  This 
new,  added  volume  can  be  obtained  by  the  expression 


(A3) 


The  new  concentration  is 


t  _  ^solute 
Vf^V/ 


(A4) 


Inserting  Eqs.  Al,  A2  and  A3  into  Eq.  A4,  we  obtain  the  expression 


(A5) 


which  reduces  to 


(A5) 
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APPENDIX  B.  Progran  ERRCHI>3.BAS 


10 

20 

30 

40 

50 

60 

70 

100 

110 

120 

130 

140 

143 

144 


REM  This  program  calculates  CHI-3  and  Gamma  for  samples  in  solution. 
REM  It  can  be  used  in  conjunction  with  Lotus  spreadsheet  CHI-3. WKl 
REM  which  uses  a  least-squares  curve-fitting  routine  to  determine 
REM  solute  CHI-3  values  from  solution  CHI-3  values  calculated 
REM  with  this  program.  B.  Kimball/  18  March  1994. 

CLS 

INPUT  "I  Of  solution  (mv)  =  ",  Isol 

INPUT  "I  of  CS2,  (mv)  =  ",  Ics 

REM  ICS  =  50000 

INPUT  "Absorbance,  =  ",  ab 

INPUT  "Concentration,  (mg/ml)  =  ",  con 

REM  PARAMETERS  FOR  CHI-3  AND  GAMMA  EVALUATION 


145  REM  Index  of  refraction  of  solution 
147  ns  =  1.44 

150  REM  Path  Length  CS-2,  (mm) 

155  DCS  =  2 

160  REM  Path  length  solution,  (mm) 

162  DS  =  2 

165  REM  CHI-3  of  Solvent,  (esu) 

167  Xsolv  =  2.4E-14 

168  REM  Uncertainty  in  Solvent  CHl-3 
170  WXsolv  =  9.7E-16 

172  REM  CHI-3  of  CS-2 

174  XCS  =  6.8B-13 

175  REM  Index  of  refraction  of  CS-2 
180  ncs  =  1.63 

182  REM  Molecular  weight  of  solute,  g/mole 
184  M  =  678 

186  REM  Avogadros  nvunber 

188  A  =  6.023E+23 

190  REM  Mass  of  sample  (mg) 

195  mg  =  3.8 

200  REM  Volume  of  Solvent  used  in  solution 
205  ml  =  5 

210  REM  PARAMETERS  FOR  CHI-3  UNCERTAINTY 

220  PWXCS  =  0 

230  WXcs  =  XCS  *  PWXCS 

240  PWIsol  =  .05 

250  WIsol  s  Isol  *  PWIsol 

260  PWIcs  =  .05 

270  WIcs  =  ICS  *  PWICS 

280  PWns  =  .01 

290  Wns  =  ns  *  PWns 

300  PWncs  =  .001 

310  Wncs  =  ncs  *  PWncs 

320  PWDS  =  .001 

330  WDS  =  DS  *  PWDS 

340  PWDCS  =  .001 

350  WDCS  =  DCS  *  PWDCS 

360  PWab  =  .0008 

370  Wab  *  ab  *  PWab 

380  PWM  =  .001 

390  WM  =  M  *  PWM 

392  PWmg  =  .01 
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394  Wmg  =  mg  *  PWmg 
396  PWml  -  .01 
398  Wml  =  ml  *  PWml 


400  REM  CALCULATE  CHI-3  AND  GAMMA 

405  alpha  =  -LOG  (10  -ab) 

410  PI  =  Isol  /  Ics 
420  P2  =  ns  /  ncs 
430  P3  =  DCS  /  DS 

440  P4  =  (alpha  /  (1  -  (EXP ( -alpha) )) )  *  EXP(alpha  /  2) 
450  Xsol  =  XCS  *  (PI  *  .5)  *  (P2  "  2)  *  P3  *  P4 
460  Xsam  =  (Xsol  -  Xsolv)  /  (con  /  1000) 

470  L=  (ns  ^2+2)  /3 

480  GAMMA  -  (Xsam  *  M)  /  (L  ^  4  *  A) 


490  REM  CALCULATE  Xsol  UNCERTAINTY 

495  REM  Sensitivity  Expressions 

500  SXcs  =  (PI  ''  .5)  *  (P2  ^  2)  *  P3  *  P4 

510  SIsol  =  XCS  *  .5  *  ((Isol  ^  .5)  ^  -1)  *  ((Ics  ^  .5)  -1)  * 

P2  2  *  P3  *  P4 

520  Sics  =  XCS  *  .5  *  ((Isol  .5)  /  (Ics  ^  1.5))  *  P2  ^  2  *  P3  *  P4 

530  Sns  =  XCS  *  (PI  ^  .5)  *  2  *  (ns  /  (ncs  2))  *  P3  *  P4 

540  Sncs  =  XCS  *  (PI  .5)  *  ns  *  2  *  (2  /  (ncs  ^3))  *  P3  *  P4 

550  SDCS  =  XCS  *  (PI  ^  .5)  *  (P2  ^  2)  *  (1  /  DS)  *  P4 

560  SDS  =  XCS  *  (PI  *  .5)  *  (P2  '•  2)  *  (DCS  /  (DS  ^  2))  *  P4 

570  Sab  =  XCS  *  (PI  ^  .5)  *  P2  ^  2  *  P3  *  ( (1  -  EXP(-alpha))  *  ((alpha  /  2)  + 

EXP(alpha  /  2))  -  alpha  *  EXP(-alpha  /  2))  /  ((1  -  EXP(-alpha))  "  2) 

580  REM  Parameter  Uncertainties 

590  UXcs  =  SXcs  *  WXcs 

600  Ulsol  =  SIsol  *  WIsol 

610  UIcs  =  Sics  *  WIcs 

620  Uns  =  Sns  *  Wns 

630  Uncs  =  Sncs  *  Wncs 

640  UDCS  =  SDCS  *  WDCS 

650  UDS  =  SDS  *  WDS 

660  Uab  =  Sab  *  Wab 

670  REM  Xsol  Uncertainty 

680  WXsol  =  (UXcs  ''  2  +  Ulsol  ''  2  *  UIcs  ''  2  +  Uns  2  +  Uncs  '•  2  + 

UDCS  ^  2  +  UDS  ^  2  +  Uab  ^  2)  ^  .5 


700 

REM 

CALCULATE  con  UNCERTAINTY 

710 

REM 

Sensitivity  Expressions 

720 

Smg 

=  1  /  ml 

730 

Sml 

=  mg  /  (ml  ^  2) 

740 

REM 

Parameter  Uncertainties 

750 

Umg 

=  Smg  *  Wmg 

760 

Uml 

=  Sml  *  Wml 

770 

REM 

con  Uncertainty 

780 

Wcon  =  (Umg  ^  2  +  Uml  ''  2)  ''  ,S 
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800  REM  CALCULATE  Xsam  UNCERTAINTIES 

810  REM  Sensitivity  Expressions 
820  SXsol  =  1000  /  con 
830  SXsolv  =  1000  /  con 

840  Scon  =  1000  *  (Xsol  -  Xsolv)  /  (con  ^  2) 

850  REM  Parameter  Uncertainties 
860  UXsol  =  SXsol  *  WXsol 
870  UXsolv  =  SXsolv  *  WXsolv 
880  Ucon  =  Scon  *  Wcon 

890  REM  Xsam  Uncertainty 

900  WXsam  =  (UXsol  ^  2  +  UXsolv  2  +  Ucon  "2)  ^  .5 


950  REM  CALCULATE  L  UNCERTAINTIES 

960  REM  Sensitivity  Expressions 
970  SLn  =  2  *  ns  /  3 

980  REM  Parameter  Uncertainties 
990  ULn  =  SLn  *  Wns 

1000  REM  L  Uncertainty 
1010  WL  =  ULn 


1020  REM  CALCULATE  GAMMA  UNCERTAINTIES 

1030  REM  Sensitivity  Expressions 
1040  SXsam  =  M  /  ((L  ^  4)  *  A) 

1050  SM  =  Xsam  /  ((L  4)  *  A) 

1060  SL  =  4  *  Xsam  *  M  /  ((L  ^  5)  *  A) 

1070  REM  Parcuneter  Uncertainties 
1080  UXsam  =  SXsam  *  WXsam 
1090  UM  =  SM  *  WM 
1100  UL  =  SL  *  WL 

1110  REM  GAMMA  Uncertainty 

1120  WGAMMA  =  (UXsam  ''2+UM^2  +  UL^2)  ^  .5 


1200  REM  PRINT  RESULTS 

1204  PRINT  "  RESULTS" 

1205  PRINT 

1210  PRINT  "CHI-3  of  solution  is";  Xsol 
1220  PRINT  "Uncertainty”;  WXsol 
1230  PRINT 

1240  PRINT  "CHI-3  of  sample  is";  Xsam 
1250  PRINT  "Uncertainty";  WXsam 
1260  PRINT 

1270  PRINT  "Gamma  of  the  solute  is”;  GAMMA 
1280  PRINT  "Uncertainty";  WGAMMA 
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PARAMETER  UNCERTAINTIES 


1310  PRINT  " 

1320  PRINT 
1330  PRINT  "  Xsol" 

1340  PRINT  "CHI-3  of  CS2  "?  UXcs 

1350  PRINT  "Voltage  reading  for  solution";  OIsol 

1360  PRINT  "Voltage  reading  for  CS2";  UIcs 

1370  PRINT  "Solution  refractive  index";  Uns 

1380  PRINT  "CS2  refractive  index";  Uncs 

1390  PRINT  "Solution  thickness";  ODS 

1400  PRINT  "CS2  thickness";  UDCS 

1410  PRINT  "Absorption";  Uab 

1415  INPUT  ;  "Hit  Return",  me 

1416  CLS 
1420  PRINT 

1430  PRINT  "  Concentration" 

1440  PRINT  "Weight";  Umg 
1450  PRINT  "Volume";  Uml 
1470  PRINT 

1480  PRINT  "  Xsam" 

1490  PRINT  "Solution  CHI-3";  UXsol 
1500  PRINT  "Solvent  CHI-3";  UXsolv 
1510  PRINT  "Concentration";  Ucon 
1520  PRINT 

1530  PRINT  "  Field  Factor  L" 

1540  PRINT  "Solution  refractive  index";  ULn 
1550  PRINT 

1560  PRINT  "  GAMMA" 

1570  PRINT  "Sample  CHI-3";  UXsam 
1580  PRINT  "Molecular  weight";  UM 
1590  PRINT  "Field  Factor  L";  UL 
1700  END 
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APPENDIX  C.  Parameter  Sensitivity  Curves 


1.7 


Fig.  C4.  Uncertainty  in  (y)  resulting  from  percent  error  in 
the  DFWM  phase^con jugate  signal  of  the  solution. 
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Pig.  C18.  Uncertainty  in  <y)  resulting  from  percent  error  in 
measurement  of  solution  weight. 
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Uncertainty  In  Gaima  C®buD 
crimes  10E-31D 


3.4 


Fig.  C20.  Uncertainty  in  (y)  resulting  from  percent  error  in 
measurement  of  solution  volume. 
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Fig.  C21.  Uncertainty  in  <y)  resulting  from  percent  error  in 
molecular  weight  of  solute* 
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COMMANDER 

U.S.  ARMY  NUCLEAR  AND 
CHEMICAL  AGENCY 
ATTN :  MONA-ZB  BLDG .2073 
7500  BACKLICK  ROAD 
SPRINGFIELD,  VA  22150-3198 


1  copy  to: 

DIRECTOR 

DEFENSE  NUCLEAR  AGENCY 
ATTN:  TDTR 

6801  TELEGRAPH  ROAD 
ALEXANDRIA,  VA  23310-3398 


1  copy  to: 

DIRECTOR 

DEFENSE  NUCLEAR  AGENCY 
ATTN:  HRP 

6801  TELEGRAPH  ROAD 
ALEXANDRIA,  VA  23310-3398 


1  copy  to: 

HQDA 

OFFICE  OF  THE  SURGEON  GENERAL 
ATTN:  DASG-HCG 

5109  LEESBURG  PIKE 
FALLS  CHURCH,  VA  22041-3258 


1  copy  to: 

COMMANDER 

U.  S.  ARMY  COMBINED  ARMS  DEFENDE  CENTER 

ATTN:  ATZL-CAD-N 

FT.  LEAVENWORTH, KS  66027-5300 


1  copy  to: 

U.S.  ARMY  AVIATION  MEDICAL  RESEARCH 

LABORATORY 

ATTN:  SGRD-VAB-CB 

FT.  RUCKER,  AL  36352-5000 


1  copy  to: 

COMMANDER 

WALTER  REED  INSTITUTE  OF  RESEARCH 
ATTN:  SGRD-UWZ 

WASHINGTON,  DC  20307-5100 


1  copy  to: 

U.S.  ARMY  RESEARCH  LABORATORY 
ATTN:  SLCHD-NW-TN 

2800  POWDER  MILL  ROAD 
ADELPHI,  MD  20783 


36 


1  copy  to 


U.S.  NAVY  CLOTHING  AND  TEXTILE 
RESEARCH  FACILITY 
ATTN;  CODE  40.1 
21  STRATHMORE  ROAD 
NATICK,  MA  01760-2490 


2  copies  to: 

DEFENSE  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 


7  copies  to; 

COMMANDER 

U.S.  ARMY  NATICK  RD  &  E  CENTER 
TECHNI CAL  L I BRARY 
ATTN;  STRNC-MIL 
NATICK,  MA  01760 


1  copy  to; 

COMMANDER 

U.S.  ARMY  NATICK  RD  &  E  CENTER 
ATTN;  STRNC-MSR 
NATICK,  MA  01760-5020 


50  copies  to: 

COMMANDER 

U.S.  ARMY  NATICK  RD  &  E  CENTER 
ATTN;  STRNC-YSD  (Mr.  Kimball) 
NATICK,  MA  01760-5020 
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